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The Diesel refrigeration cycle using an ideal quantum gas as the working substance is called quantum 
Diesel refrigeration cycle, which is different from Carnot, Ericsson, Brayton, Otto and Stirling refrigeration 
cycles. For ideal quantum gases, a corrected equation of state, which considers the quantum behavior of 
gas particles, is used instead of the classical one. The purpose of this paper is to investigate the effect of 
quantum gas as the working substance on the performance of a quantum Diesel refrigeration cycle. It is 
found that coefficients of performance of the cycle are not affected by the quantum degeneracy of the 
working substance, which is the same as that of the classical Diesel refrigeration cycle. However, the 
refrigeration load is different from those of the classical Diesel refrigeration cycle. Lastly, the influence 
of the quantum degeneracy on the performance characteristics of the quantum Diesel refrigeration cycle 
operated in different temperature regions is discussed. 

© 2008 Elsevier Ltd. All rights reserved. 


1. Introduction 

Quantum theory, a field opened by Planck and other pioneers in 
20th century has intrigued scientists. Quantum thermodynamic 
cycles, such as quantum engine [1-4], quantum refrigerator [5], 
quantum heat pump [6], quantum amplifier [7], quantum after¬ 
burner [8], and so on [9-12], have become one of the interesting 
research subjects for people working in thermodynamics and sta¬ 
tistical physics. These quantum analyses provide the application 
foundation of the equilibrium or non-equilibrium statistical 
mechanics to the practical engineering cycles. Several authors have 
intensively investigated the performance characteristics of quan¬ 
tum thermodynamic cycles working with particles in a potential 
[1], uncoupled spins [2], harmonic oscillators [3], two level sys¬ 
tems [5], multi-level systems [4] or quantum optics systems [8- 
10], etc. 

The Diesel cycle is one of the important cycle models in the re¬ 
search and manufacture. Based on the equation of state of an ideal 
classical gas, many novel results of the Diesel cycle are obtained 
[13-19]. Recently, the thermodynamic performance parameters 
of the Carnot [20], Stirling [21-23], Ericsson [24], Brayton [25,26] 
and Otto [27-29] cycles working with an ideal quantum gas have 
been investigated by a number of researchers. The performance 
characteristics of these cycles are analyzed under the condition 
of the quantum degeneracy. However, most of the previous works 
have concentrated on the performance analysis of the above these 
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cycles. The performance analysis of the Diesel refrigeration cycle 
working with an ideal quantum gases has been rarely considered. 

In this paper, the quantum models of Diesel refrigeration cycle 
working with ideal Fermi ( 3 He) and Bose ( 4 He) type quantum gas 
are established, which are called Bose or Fermi Diesel refrigeration 
cycle, respectively. Based on the corrected equation of state of an 
ideal quantum gas, the general expressions of the coefficient of 
performance and refrigeration load of the cycle are derived. The 
influence of the quantum degeneracy on the performance para¬ 
meters of the cycle is discussed. The results of the work obtained 
here will be compared with other quantum gas refrigeration cycles. 


2. A gas Diesel refrigeration cycle 

A Diesel refrigeration-cycle using an ideal quantum gas as the 
working substance is composed of one isobaric, one isochoric and 
two isentropic processes. According to the expression of the entro¬ 
py of an ideal gas, one can easily give the T-S diagram of a Diesel 
refrigeration cycle, as shown in Fig. 1. Process 1-2 is an isentropic 
compression. Processes 2-3 and 4-1 are isobaric heat rejection and 
isochoric heat addition, respectively. Process 3-4 is an isentropic 
expansion. Qh and Ql are the amounts of heat exchange between 
the working substance and the heat reservoirs at temperatures 
T h and T l during isobaric 2-3 and isochoric 4-1 processes; where 
p H is the highest pressure of the working substance in the cycle, 
v H is the maximum average volume that a particle occupies. T H 
and T l are the temperatures of the hot reservoir and cold reservoir. 
The temperatures of state points 1, 2, 3, and 4 are represented by 
Ti, T 2 , T 3 and T 4 , respectively. 
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Nomenclature 



Cp 

heat capacity at constant pressure (J/K) 

S 

entropy (J/K) 

Cv 

heat capacity at constant volume (J/I<) 

T 

absolute temperature (K) 

F(z) 

correction function 

T h 

temperature of heat reservoir (I<) 

g 

weight factor 

T l 

temperature of cold reservoir (K) 

h 

Planck’s constant (Js) 

u 

internal energy (J) 

k 

Boltzmann constant (J/K) 

V 

volume (m 3 ) 

m 

rest mass of a particle (kg) 

w 

work input (J) 

N 

total number of particles 

z 

fugacity of the gas 

N 0 

particle number in the ground state 



n 

number density (m -3 ) 

Greek 

symbols 

P 

pressure (Pa) 

s 

coefficient of performance 

Ph 

maximum pressure (Pa) 

X 

mean thermal wavelength (m) 

Qh 

amount of heat exchange between the working sub- 

P 

chemical potential of the gas (J) 


stance and the heat reservoir during the high con- 

T 

temperature ratio of two heat reservoirs 


stant-pressure process (J) 

k 

Boltzmann’s constant 

Ql 

amount of heat exchange between the working sub- 

m 

Gamma function 


stance and the heat reservoir during the high constant- 

cry 

Riemann-Zeta function 


volume process (J) 



Rql 

relative refrigeration load 




3. Thermodynamic properties of an ideal quantum gas 


Based on quantum statistics, the expressions for the pressure 
and number density of an ideal Fermi gases are given as follows 
[30] 

P = s«W) (1) 


= N = 1 gf 3/2 (z) 

~v~v~ y 



where/i(z) = r ], J 0 °° dx is called the Fermi function and g is a 
weight factor that arises from the “internal structure” of the parti¬ 
cles (such as spin), k is the Boltzmann’s constant, T is the gas tem¬ 
perature, N is the total number of the particles, and V is the volume 
of the gas, v is the occupied volume of the particle, 
X = h(2nmkT) 1/2 is the mean thermal wavelength of the particles, 
h is Planck’s constant, m is the rest mass of the particles, //(z) is 
called the Fermi function, z = exp [p/kT] is the fugacity of the gas, 
r(l) is the Gamma function and /i is the chemical potential of the 
gas. 

For an ideal Bose gas the expressions for the pressure and the 
number density are given as follows [30] 


kTg s/2 (z) 

A 3 




Fig. 1. Temperature-entropy diagram of a Diesel refrigeration cycle. 


n 


N - No 
V 


1 = £3/2(2) 
v 


X 1 


(4) 


where N 0 is the number of the particles in ground state, while 


OO x l-1 

r(l) JO 


dx is called the Bose function. 


Si(z) = no Jo 

The equation of state for an ideal quantum gas may be written 
as [24] 

P = nkTF(z) (5) 

where F(z) is called the correction factor which is defined for a Fer¬ 
mi gas as 

fs/2 (z) 


F(z) = 


F(z) 


f3/2 (Z) 

The correction factor for a Bose gas is 

g 5/2 (z) /n 3/2 


( 6 ) 


gi/2 (z) \Tc 


T x 3/2 


- 1 


0(T - T c ) + 1 


0(T - T c ) = 


1 (T>T C ) 


(7) 


0 (T < T c ) 

where T c is the Bose-Einstein condensation temperature defined as 


W 


2nmk 


N 


_ 2 

3 


w 


VC(3/2)J 2nmF V ^ 3 ^ 


( 8 ) 


where £(3/2) = 2.612 . 

The internal energy and the specific heat at constant volume of 
an ideal quantum gas are given by 


U = -NkTF(z ) 


C 


v 


'0LT 

0T 


v 


Cp = Cy + T 


1mJL[TF(T, v)) 




/0P 

\vrj v \dT yp 


(9) 


( 10 ) 


( 11 ) 


Substituting Eqs. (5) and (10) into Eqs. (11), one can obtain 
Cp = 2 Nk d7 [TF(T ' P)1 


(12) 
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where the correction factor F(T, v) is a function of temperature T and 
volume v, the correction factor, F(T, p ) is a function of temperature T 
and pressures p. Using Eqs. (11) and (12), one can calculate the 
amounts of heat exchange in isochoric and isobaric processes. They 
are given by 

a,; = J C V (T , v)dT = | Nk[T f F(T f , v) - T,F(T,, v)] (13) 

and 

Q£ = f C P (T, p)dT = | Nk[T f F(T f , p) - T t F(J u p)] (14) 


Using Eqs. (16), (18), (21), and (22), yields 




\NkT H F(T H ,p H ) 




Using Eqs. (19) and (23), one can obtain 




*NkT H F(T H , p H ) 


V 2 

V 3 


1 


5. Expressions for several important parameters 


(25) 


(26) 


where the indices i and /refer to the initial and final states, respec¬ 
tively. We can analyze the Diesel refrigeration cycle working with 
an ideal quantum gas with the help of Eqs. (13) and (14). 

4. The effect of quantum degeneracy 

Since the entropy s of the system and the correction factors F(z) 
only depend on z, the expressions in the isentropic processes 1-2 
and 3-4 can be obtained 


The net work is easily obtained by Eqs. (25) and (26) as follows: 

W = Q h - Ql 



NkT H F(T Hl p H ) 


1 


3 

5 




(27) 


Using Eqs. (26) and (27), the coefficient of performance (COP) is 
written as 


F(T Ll Vh) — F(T 2 , Ph) 


F(Th, Ph) = F(T 4 , Vh) 


(15) 

(16) 



1 


Qh i 
Ql 1 


1 


5 7 , 3 V2~ V 3 
3 5 5 

^ 2-^3 


1 


(28) 


Because of this property of the correction factors, quantum and 
classical ideal-gases have the same isentropic relations between 
temperature, pressure and specific volume. Consequently, for the 
isentropic processes 1-2 and 3-4, one can obtain the following 
equations: 


V H _ fpA 5 _ f[l ) 2 

t' 2 ~[pj ~\Tl) 



Eq. (28) is the same as that of the classical Diesel refrigeration 
cycle s c = —j—*-. It is understood that the COP value of the Diesel 

IT O V') — Vo 


refrigeration cycle is not affected by the quantum degeneracy of 
the refrigerant. In addition, in order to compare the performance 
of a quantum Diesel refrigeration cycle with that of a classical Die¬ 
sel refrigeration cycle, we introduce the relative refrigeration load 
as follows 


v 2 3 \pj \T 4 J 


(18) % = § = F(T„, p H ) 



Using Eq. (5), for the isobaric process 2-3 and the isochoric 
process 4-1, one can obtain 


^2 _ Vi F(T h , p H ) 
T h v 3 F(T 2 ,P h ) 


T 4 _ P 4 F(Tl, Vh) 

Fl Pi F(T 4 , Vh) 

Using Eqs. (17) and (19), yields 

n = (Tl\ (IA (vA 2 ' 3 V 2 FJTh, p H ) 
Th ViJVh) \v h ) V 3 F(T 2 , Ph ) 




The amounts of heat exchange, Ql during isochoric 4-1 process 
and Qh during isobaric 2-3 process, are obtained 

Ql = f' C - dT = f L C ” dJ = v H ) - T 4 F(T 4 , v H )i 

JT 4 JT a Z 

( 22 ) 


or 


Rql = % = F(T l , Vh) 


where Q[ = \NkT H 


(30) 


V 2 

Vh 


V 2 

v 3 


”3 


21 


or Q( = | NkT L 




53 


IS 


the refrigeration load of a classical Diesel refrigeration cycle. 


6. Several interesting cases 

It is clearly seen from Eqs. (29) and (30) that, for a quantum Die¬ 
sel refrigeration cycle, the relative refrigeration load is dependent 
on the temperature, volume, pressure and other parameters. It is 
obvious that there are different performance characteristics for 
the quantum Diesel refrigeration cycle operated in the different 
temperature regions. Some special cases are discussed. 

6.1. Under the condition of strong gas degeneracy 





^Nk[T 2 F(T 2 ,p H )-T H F(T H ,p H )} 



where T H = T 3 , T L = T u v { = v 4 = v H and p 2 = p 3 = Ph (see Fig. 1 ). Using 
Eqs. (16), (20), and (22), yields 


Q_ l = J t C v dT = ^NkT L F(T L , v„) 


1 - 


^2 

V 3 



In this case, the working gas is always completely a degenerate 
quantum gas throughout the cycle. The correction factor F(z) of an 
ideal Fermi gas can be obtained in terms of v as [22,23] 


F(T, v ) 


2 T f (v) k 2 T 

5 T + 6T f (^) 


(31) 


2 3 ~ 

where T F (v) =Cv~ 5 is the Fermi temperature, C = (3h 2/ ^ 7r)3 is a con¬ 
stant. By using Eq. (31), the relative refrigeration load Eq (30) can 
be rewritten as follows 
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Rql — F(Tl, Vh) 


2C V \\ 
5 T l 


n 2 T L 
6 Cv^ 



Similarly, in this case, the correction factor F(z) for an ideal Bose 
gas can be obtained in terms of v as 


F(T, v) = f 



f2nmk\ 2 


h d 


1 


vT 2 



Where £(§) = 1.341. By using Eq. (33), the relative refrigeration 
loads Eq. (30) can be rewritten as follows 


Rql = F(Tl , Vh) = C 



2nmk\ 
h 2 ) 




6.2. Under the condition of weak gas degeneracy 


In this case, the working gas is in the high-temperature and 
low-density condition throughout the cycle. The correction factor 
F(z) of an ideal Fermi gas can be obtained in terms of v as [22,23] 

F(T,v) = 1+3E (35) 

vV 

3 

where G = —-—j is a constant. By using Eq. (35), the relative refrig- 

32(nmk)2 

eration loads Eq. (30) can be rewritten as follows 


Rql=F(T l , v h ) = 1+-^ t (36) 

VhTI 

Similarly, in this case, the correction factor F(z) for an ideal Bose 
gas can be obtained in terms of v as 

F(T, v) = 1 - (37) 

vT 2 

3 

where G = —is a constant. By using Eq. (37), the relative refrig- 

\6(Kmk)2 

eration load Eq. (30) can be rewritten as follows 


Rql = F(T l , v„) = 1-(38) 

When the temperature of the gas is high enough and its density 
is low enough, the fugacity of an ideal quantum gas z is much smal¬ 
ler than unity. In such a case, F(z) = 1 and an ideal quantum gas be¬ 
comes an ideal classical gas. The relative refrigeration load can be 
expressed as follows 

Rql = 1 (39) 

The results are just those of a classical Diesel refrigeration cycle. 
This shows clearly that, at high temperatures, the quantum Diesel 
refrigeration cycle becomes the classical cycle. 


7. Discussion 

In Fig. 2, it is shown how the relative refrigeration load R F QL de¬ 
pends on the temperature T L and the maximum occupied volume 
% of the particle in the case of Fermi gas. It shows clearly that un¬ 
der the condition of strong gas degeneracy, the refrigeration load 
Fql is higher than that of weak gas degeneracy. For high values of 
temperature T L , the relative refrigeration load R F QL goes to unity. 
This is the expected behavior since correction factors goes to unity 
for high values of temperature T L and the refrigerant becomes a 
classical ideal gas. 

In Fig. 3, it is shown how the relative refrigeration load R B QL 
depends on the temperature T L and the maximum occupied 
volume v H of the particle in the case of Bose gas. It shows clearly 



0 



Fig. 2. The three-dimensional graph of the relative refrigeration load R ^ varying 
with the cold reservoir T L and the maximum occupied volume v H of the particle in 
the case of Fermi gas. 


the refrigeration load R B QL under the condition of strong gas degen¬ 
eracy is lower than that of weak gas degeneracy. For high values of 
temperature F L , the relative refrigeration load R B aL goes to unity. 
This is also the expected behavior since correction factors goes to 
unity for high values of temperature T L and the refrigerant be¬ 
comes a classical ideal gas. 

The relative refrigeration load and R^ varying with the tem¬ 
perature T l and the maximum occupied volume v H of the particle 
are shown in Figs. 2 and 3. For low values of J L and v H} the relative 
refrigeration load R F aL and R B QL are different from the unity since 
quantum degeneracy becomes important. It is seen that R F aL is al¬ 
ways greater than unity and R^ is always lower than unity. On 
the other hand, the relative refrigeration load and R^ go to 



1 

0.75 

Rql 0.5 

0.25 

0 

lx 


Fig. 3. The three-dimensional graph of the relative refrigeration load Rq L varying 
with the cold reservoir T L and the maximum occupied volume v H of the particle in 
the case of Bose gas. 
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unity for high values of T L and v H . It also can be seen that the rel¬ 
ative refrigeration load of Diesel refrigeration cycle can be ex¬ 
pressed as follows 

4 > 4 = 1 > 4 ( 4 °) 

Moreover, unlike Bose-Brayton refrigeration cycle, the mini¬ 
mum value of the cold reservoir T L for the Bose Diesel refrigeration 
cycle is not restricted by the Bose Einstein condensation tempera¬ 
ture T c . 

8. Conclusion 

In this paper, a thermodynamic model of the quantum Diesel 
refrigeration cycle has been established. The influence of the quan¬ 
tum degeneracy of the working fluid on the performance of a Die¬ 
sel refrigeration cycle is analyzed. Under the quantum degeneracy 
conditions, it is shown that the refrigeration load of the Fermi Die¬ 
sel cycle is always greater than that of the classical Diesel cycle. On 
the contrary, the refrigeration load of the Bose Diesel cycle is al¬ 
ways lower than that of the classical one. Thus, in a Diesel refriger¬ 
ation cycle using a Fermi-type refrigerant, more heat is absorbed 
per cycle from the cold reservoir by consuming more work per cy¬ 
cle. On the other hand, in a Diesel refrigeration cycle using a Bose- 
type refrigerant, less heat is absorbed per cycle from the cold res¬ 
ervoir by consuming less work per cycle. It can be said that the Fer¬ 
mi type of ideal gas is always advantageous for a Diesel 
refrigeration cycle since the same refrigeration load is obtained 
with a lower number of cycles. On the contrary, the Bose type of 
ideal gas is always disadvantageous because it causes an increment 
in the number of cycles for a given refrigeration load. Moreover, 
unlike Bose-Brayton refrigeration cycle, the minimum value of 
the cold reservoir T L for the Bose Diesel refrigeration cycle is not re¬ 
stricted by the Bose Einstein condensation temperature T c . 
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